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I. INTRODUCTION 



Establishing the vahdity of various cosmological scenarios requires understanding the 
dynamical evolution of the cosmos in time, e. g., the time evolution of the inflation-driving 
field |jl|] before, during and after inflation, and the details of symmetry changing phase 
transitions that may have given rise to cosmic strings 0. The conventional formulation 
of quantum field theory in terms of causal Green's functions in the Heisenberg picture is 
not especially suited to time- dependent problems that make use of an initial condition for 
specific solution. Green's functions contain all information needed for determining transition 
rates, 5'-matrix elements, etc., of systems in equilibrium where initial data are superfluous. 
However, following a system's time evolution from a definite initial configuration is more 
efficiently accomplished in a Schrodinger picture description, where the initial data consist 
of specifying a pure or mixed state. 

For bosonic fields, the field-theoretic functional Schrodinger picture is a generalization 
from ordinary quantum mechanics to the infinite number of degrees of freedom that con- 
stitute a field. Therefore, it allows the use of the mathematical/physical intuition acquired 
in quantum mechanics to analyze field-theoretic problems. Notwithstanding, the functional 
Schrodinger picture is not as widely used in actual calculations as the Green's function 
method since renormalization is more easily carried out in the latter framework. However, 
it has been established renormalizability for of the Schrodinger picture for both static Q 
and time-dependent p cases. 

In the Schrodinger picture, a pure state is given by a single wave functional of a 
c-number field 0(x) at fixed time, while a pure state is described by a functional density 
matrix p(0i,02) = J2nPn'^n{4>i)'^n{4>2), whcrc {"^n} is a complete set of wave functional 
and Pn is the probability that the system is in the state n |0. The time evolution of 
pure and mixed is governed by the time-dependent Schrodinger equation and the Liouville- 
von Neumann equation, respectively. In general these equations cannot be solved, except 
for systems whose Hamiltonian is quadratic, and consequently approximation methods are 
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needed. Variational approximations were developed in Refs. which lead to tractable 

equations in the cases of Gaussian trial states. The Gaussian approximation leads to self- 
consistent equations that, unlike perturbation theory, reflect some of the non-linearities of 
the full quantum theory. 

In various inflationary models, the dynamics of the early universe is dictated by the time 
evolution of an scalar field, usually called inflaton. The equation of motion used in most of 
the analyses is of the form 

^ + 3-0 + Kff = , (1) 

a 

where Ves is the static, finite-temperature effective potential. However, the static effective 
potential does not properly take into account effects of non-equilibrium dynamics. Moreover, 
semiclassical or linearized approximations have been frequently made, and the full 



non-linearity of an interacting theory is lost. It is conceivable that a more complete analysis 
may change the picture of the early universe drawn from these approximations |]TT] , p]2| . 

In this work we study a self-interacting scalar quantum field model in a flat Robertson- 
Walker space-time in the functional Schrodinger picture. Our purpose is to obtain a set of 
quantum dynamical equations describing the evolution of a scalar field as well as the scale 
factor. In order to obtain a small set of tractable equations, which capture some of the non- 
linearities of the full problem, we employ a variational method using a Gaussian trial state 
whose kernels depend on just a few parameters. We analyze the renormalization of the equa- 
tions of motion for the variational parameters to assure the consistency of them. Moreover, 
we also study the renormalization of the expectation value of the energy-momentum tensor, 
which is the source for the semi-classical Einstein equation, showing that it is finite once 
we take into account the usual mass and coupling constant renormalizations. Therefore, we 
obtain a set of equations that can be used to address various cosmological questions like the 
dynamics of chaotic inflation. 

This article is organized as follows: In Sec. II we present the functional Schrodinger 
picture for a scalar field and obtain the Gaussian approximation for the ground state of the 



model, which will be used to define our simplified Gaussian approximation. Section III 
contains the Ansatz that we used in our calculations as well as the study of the renormal- 
ization of the equations of motion for the variational parameters. We present the analyses 
of the renormalization of the energy-momentum tensor in Sec. IV and, in Sec. V we draw 
our conclusions. 



II. FUNCTIONAL SCHRODINGER PICTURE 
A. Generalities 

Here we present the basic facts about the functional Schrodinger picture, and the in- 
terested reader can find exphcit examples and learn more about it in Refs. [^|6|,p!T|. In the 
field-theoretic Schrodinger picture, pure states are described by wave functionals of a 
c-number 0(x) at a fixed time. The inner product is defined by functional integration 

while operators are represented by functional kernels. 

0|^) ^ j D(j)' O(0,0')^(0') (3) 

We adopt a diagonal kernel ^'(x) =^ (j){x.)5{(j) — (j)') for the canonical field operator $(x) at 
fixed time, while the canonical commutation relations determine the canonical momentum 
kernel to be n(x) =^ {l/i)[6/6(f){x)]6{(f)—(j)'). Hence $ acts by multiplication on functionals of 
and H acts by functional differentiation. In this way, the action of any operator constructed 
from H and $ is 



C(n, $)|^) ^ O 



1_S_ 

i 5(b 



^(0) • (4) 



The fundamental equation for initial value problems is the time- dependent Schrodinger 
equation for the time-dependent state functional ^1^(0; t). This equation takes definite form, 
once a Hamiltonian operator H{Il, $) is specified: 



i 64> 



vl/(0;t). 



(5) 



The initial value problem is completely defined once we also supply the initial wave func- 
tional. 



B. Dirac's Variational Principle 



The time-dependent Schrodinger equation (||) cannot be directly integrated, unless the 
system is described by a quadratic Hamiltonian. Therefore, we shall employ a variational 
approximation to study non-linear (interacting) systems. Applications of variational prin- 
ciples with restricted variational Ansatz, in the Rayleigh-Ritz manner, result in tractable 
self-consistent dynamical equations for the parameters used in the Ansatz, which still retain 
some of the non-linearity of the complete problem. 

For pure states, the time-dependent Schrodinger equation can be obtained through 
Dirac's variational principle [0. First we define the effective action F as the time inte- 
gral of the diagonal matrix element of id/dt — H: 

d 



dt m 



1—-H 
dt 



(6) 



and then we demand that F be stationary against arbitrary variations of |\E') and (^|, 
imposing appropriate boundary conditions. 

It is possible to implement Dirac's principle in two steps in such way that F can be 
associated to the functional generator of the one-particle irreducible Green's functions with 
arbitrary energy and momentum |^. In order to do so, we consider the time integral of an 
off-diagonal matrix element 



F 





d 


j dt 





(7) 



subject to the constraint that the matrix element of the field $(x) is held fixed at a prescribed 
function (y9(x, t). 
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(*_|$(x)|*+) = (^(x,i) (8) 
(*-l*+) = l (9) 

These constraints are supplemented by the boundary condition that these states tend to 
the ground state of H iox t ^ ±00 The physical theory is recovered when we remove the 
constraints by solving 

5V 



. (10) 



C. Gaussian Ansatz 

Dirac's variational principle can be used to obtain approximations provided we restrict 
the variation of the trial wave functional ^' to a subspace of the full Hilbert space. In this 
work, we shall use Gaussian trial states, whose most general expression is 

t) = N{t) exp |i £ 7r(x, t) [(/)(x) - (^(x, t)] | 

xexp|-£^[0(x)-(^(x,t)] [^n-i(x,y,t)-iE(x,y,t)] [0(y) - <^(y, , (11) 

where the variational parameters are ip, tt, fl, and E and we abbreviated the integral in 
d spatial dimensions as = / (i'^x. The physical meaning of the parameters of this wave 
functional can be inferred from linear and bilinear averages. The linear averages are given 

by 

($(x)) = v.(x,0, (12) 
(n(x)) = 7r(x,i), (13) 

while bilinear averages are 

($(x)$(y)) = ifi^My) + l^(x, y, t) , (14) 
(n(x)n(y)) = 7r(x)7r(y) + ^Q-^(x, y, i) + 4 (EQE) (x, y, t) , (15) 
($(x)n(y)) = ^5(x - y) + 2 (l^E) (x, y, t) , (16) 



where we have used the matrix notation (C/C) (x, y) = J^, C(x, z)/C(z, y). Moreover, from 
the average value of the operator i{d/dt) appearing in the effective action we find that 
the imaginary part of the covariance function (S) plays the role of a canonical momentum 
conjugate to the real part Q. 

= Jj{^,t)ip{^,t) + S(x,y,t)ll(y,x,t) (17) 



D. Gaussian Vacuum in Flat Space-Time 



In order to gain some intuition and also to motivate the simplification of the Gaussian 
Ansatz, that we shall use in this work, it is interesting to obtain the vacuum state, in the 
Gaussian approximation, for the X^'^ model in Minkowski space-time. The dynamics of this 
system is governed by the Hamiltonian 



(18) 



where the potential function t/($) is 



^ ^ 2 4! 



(19) 



Substituting the wave functional (|TT]) into the Dirac's variational principle leads to 



r(y?,7r,fi,S) = J dtj 



+ n 
1 



(x, X, t) - -n'\^, X, t) - 2 (sr^s) (x, X, t) 



-(-V^fi(x,y,t) U=y +f/(2)(v9)r](x,x,t) 
- ^t/(^)(v.)fi(x,x,t)n(x,x,t)| , 



(20) 



where U^"'^ = dJ^U /dip'^. The terms in the first square bracket are the classical action, while 
terms in the second (last) square bracket are formally 0{h) {0{h'^)) corrections. In fact, 
the expression (^0|) contains all powers in h since the kernel VL must satisfy a self-consistent 
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equation, as shown below. By varying ( pOD with respect the parameters in our Ansatz we 
obtain four variational equations 

6r 



5V3(x,t) 

6n{-x, t) 

(5S(x,y,t) 
5V 

5f2(x,y,t) 



=^ 7r(x, t) = 0(x, t) 



^ TT 

= 







^](x,y,^) = 2[(r]S)(x,y,t) + (Sr])(x,y,t)] 

1 



S(x, y, t) = ^l]-2(x, y, t) - 2S2(x, y, t) 

+ + it/(^)(^)fi(x,x,t)] <5(x - y)} . 



(21) 
(22) 
(23) 

(24) 



Translation invariance implies that (p is homogeneous and that the kernels can be ex- 
pressed as a Fourier transformation (FT) 



fi(x,y,t)= / e^'^ (^-y)fi(k,t), 



(25) 



where the momentum-space integral (/ d'^'k/ {2'kY) is denoted by /j^.. Moreover, the kernels 
for the vacuum state are time-independent and the above equations of motion reduce to 





= 0, 




= 0, 


n(k)S(k) 


= 0, 


Jm) 


= 0, 



(26) 
(27) 
(28) 
(29) 



whose solution is ^ = S(k) = and 
with satisfying the gap equation 



2Vk2 + m2 ' 



(30) 



(31) 



Since this last equation is a self-consistent one for m^, some of the non-linearities of the 
complete problem are retained by the Gaussian approximation. 
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E. Renormalization of the Effective Potential 



At this point it is instructive to study the renormahzation of the Gaussian effective 
potential, which can be obtained from the effective action through 

r(y,,7r,S,fi)|^,^,.^ = -\/eff(^,7r,S,fi) / . (32) 

Therefore, from Eq. ( PUj) we can see that the Gaussian effective potential is 

+ In-'M + 2 (mE) (x,x) - ^V^ fi(x,y)|,^y 



2 

+ ^ ( /i^ + ^y?^ ) ^^(x, x) + ^fi(x, x)fi(x, x) . (33) 



The effective potential Ves{f) is obtained from V^sif, tt, S, Q) by minimizing with respect 
to the parameters vr, Q, and S. This procedure leads to vr = S(k) = and Eqs. (^) and 
(^), resulting in 



+ 



+ -/ . ' (34) 

In the limit = 3, the above integrals are clearly divergent, due the short distance 
behavior of fi, and we must renormalize V^g. We regularized these integrals by using dimen- 



sional regularization |n| on the spatial dimension (i, obtaining that the regularized effective 
potential is 



A 1 fm'^V'^'^^^'^ f3-d' 

^^'^""^ = T^' + 4! + (4vr)(^+i)/2(l-rf) [a^ ) ^ 
2m^ - / o A 
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where A is an arbitrary mass scale. Notice that the term XQQ gives rise to a double pole 
in (i = 3, while the other divergences are single poles. Since the Gaussian approximation 
is very similar to the large- approximation, the effective potential becomes finite by the 



renormalization prescription |15 

A Xr 
1 1 



1 



A Xr (47r)('^+i)/2 (d-l){3~d) ■ 
This last relation implies that 

167r2(3 - d) 



(36) 
(37) 



limA = -1Qtt\3 -d){l 



A; 



+ 0[{3-dy] \ for Xr^O. 



(38) 



When the above renormalization prescription is substituted into Eq. (p5D, it leads to a 
finite expression for the effective potential 



m 



m m 

7 + TT— ^In 



2 647r2 ' ' 647r2 \^47rA2 / 

in the limit c? = 3, where we used that 



A, 



2A 



R 



(39) 



m 
A2 



2\ ('i-3)/2 



r 



'3 - d' 



{3-d) 



+ 7 - In 



m 
A2 



0{d-3) 



(40) 



with 7 being the Euler constant. At this point we choose the scale A to be 

A^ 

which leads to 



4:71 



(41) 



m 



m 



2 "'^ ^ 647r2 





1' 








~ 2_ 



2A 



R 



2Xr 



(42) 



which is the standard result. 

The renormalized expression for m'^ can be obtained either by substituting the renor- 
malization prescription into Eq. (pT]), or by minimizing the renormalized V^g with respect 



to m . 



dm? 



^ m 



A, 



A 



R 



2 , '^R 2 , 

^« + + 32vr2 



m 



m In 



(43) 
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III. SCALAR FIELD EQUATIONS OF MOTION IN ROBERTSON- WALKER 

SPACE-TIME 



In this section we obtain the renormahzed equations of motion for a self-interacting scalar 
field using a simplified Gaussian approximation. We consider a flat Robertson- Walker space- 
time with the line element 

ds^ = df - a\t)d^\ (44) 

where a{t) is the scale factor. We assume minimal coupling between gravity and the scalar 
field and that the scalar field dynamics is governed by the Lagrange density 



C = a' 



(45) 



where the potential U is given by Eq. ([T9|). Although we are mainly interested in physical 
space-time dimensionality, d = 3, we consider the theory in d spatial dimensions in order to 
regularize the theory in later discussions. 
The canonical momentum 11 is defined by 

n = ^ = a'^^ , (46) 

and the Hamiltonian density of the system is 

Ti: = a'^ |i [a-^^U^ + a-2(V$)2] + t/(<I>)| . (47) 

The use of the general Gaussian Ansatz ( pT]) leads to coupled integro-differential equa- 
tions for the Fourier modes of the kernels Q and S Therefore, we must solve an infinity 



(large) number of coupled equations either analytically or numerically when we apply this 
approximation to study a physical problem. In order to reduce the number of free parame- 
ters and equations, we introduce a simplified Gaussian Ansatz, which is obtained by fixing 
the functional dependence on k of the kernels Q and S. In this work, we assume that these 
kernels, appearing in the Gaussian trial state ([TT]), have a form similar to their static ones, 
that is 

12 



fl{k,t) 



and 



S(k,t) 



2Jk^ + a^a{t) 



(48) 



(49) 



8[P + a2a(t)]" ' 

where a and (3 are the variational parameters and n is conveniently chosen to control the 
infinities in the approximation. Clearly, the above Ansatz recovers the vacuum solution 
in the static limit, provided we take a = 1, (3 = and a = m?. 



The effective action (y) evaluated in the Gaussian trial state ([TT|) with the kernels (^ 
and (|49|) is 



dt 



7ry9 — a 



1 /;2 \ 

-a vr H 09 H — -09 

2 2^4!^ 



k (F + a^aY+^l"^ 



+ —a^(3{2Ha + a) f 
32 Jk 



1 



k (A;2 + a2a)"+3/2 



64 ^ A 



k (A;2 + a2a)2n+i/2 



""^ / \//c2 + a^a + ^aa / , 
2 ik 4 Ay/WT 



1/2 ^2 

- -a /i + -v? 



X 



.2-d 



(50) 



32" A,k' VF + o^c^ + J 
where H = a/a is the Hubble constant. In the limit d = 3 infinities appear in F, due to 
the short distance behavior of the kernels Vl and S. However, choosing > 1 we do not 
introduce further infinities besides the one appearing in V^g ||16|| . 

Using dimensional regularization, we can separate the divergent and finite parts of F as 



dt 



a \ {'^-3)/2 
(47r)('^+i)/2(l -rf) VA2, 



.2 , '^,„2 



'l-d) 2 



A 

+ 77 



a 



a ■^id--i)/2^ 1^ 3 -rf ^ 



+ terms finite at ci = 3 



(51) 



2 (47r)('^+i)/2(l -d) VA2; ~ V 2 

The divergent part of F is similar to the divergences encountered in V^fj, see Eq. (^). This 
allow us to conclude that F can be made finite in the limit d = ?> using the renormalization 
prescription (|36|)-(|37D. It is straightforward to verify that F becomes finite in the limit c? = 3 
by this renormalization prescription: 
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dt 



TT(p —TT 



o o/ 



64 



a 



a 



2^ 647r2 



2 y 2A/J 2A/j_ 
where we chose the scale A as in Eq. (^1]) and we defined the integrals 



(52) 



1 



:{a a 



r(j - d/2) 



(53) 



/k (F + a^ay 2"'7r^/2 r(j) 

At this point it is interesting to compare our results with the ones for general Gaussian 
Ansatz as shown in Ref. |^. In both approximations the dynamical equations of motion 
become finite by the vacuum sector renormalization prescription. Moreover, the large k 
behavior of the kernels in these two approximations are similar: in the general Gaussian 
Ansatz it is required that Q ~ 0{k~^) and Q 0{k^^), while in our approximation fl 
exhibits the same high energy behavior, by construction, and n > 1, that means that the 
asymptotic k dependence of S is similar in both approximations. 

By varying the renormalized effective action F with respect to the parameters ip, vr, a, 
and (3, we obtain four coupled variational equations 

6r 



6ip 
6T 

ST 

6(3 



(5r _ ^ 

6a 









TT 



(p = a vr , 

-3^-^2n+l/2 „ 

a = a (3— 6 



n 



'n+3/2 



n — 1 



-Ha 



(54) 

(55) 
(56) 



/3 = 



,(4n + 1) (3n- 1) 2 



71+3/2 



32a^ 



V? a . 

— H t; In 

2 327r2 



2n- 1 



a 



-^2n+3/2 



{a - fil) 



Mr, 



X 



R 



n — 1 



(57) 



Due to the choice of our trial state, the equations of motion for the parameters cp and vr are 
free field ones with a time-dependent mass a, whose dynamics is dictated by the last two 



equations [[17|] . 
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IV. RENORMALIZING THE ENERGY-MOMENTUM TENSOR 



In order to write the semi-classical Einstein equation we must study the renormalization 
of the expectation value of the energy-momentum tensor in our trial state. The energy- 
momentum tensor for the scalar field described by the Lagrange density (|45|) is 

2 51 



T — 



2 ^ 2 4! 



(58) 



where / is the action, G^y = R^y — \g^uR is the Einstein tensor and the notation ; ix 
denotes the covariant derivative with respect to the space-time index /z. In the functional 
Schrodinger picture, we express the energy-momentum tensor operator in terms of the field 
operator $(x) and its canonical momentum n(x) and evaluate the expectation value in a 
given state. 

In the flat Robert son- Walker metric (H), the expectation value of T^^ in the transla- 
tionally invariant Gaussian state ( ]TT| ) has the form 



(Too) 



'2d 



-2d 



1 

+ 2 
A 



x=y 



— Jl(x, X, i)f)(x, X, t) 
8 



(59) 



-2d 



-TT 



_ ^ A ^ ^ ^ 2a-'^(Er]S)(x, X, t) 



-2d 



4! 



- ( 1 - ^) a'^Wl + /i' + ^if^ + ^^](x, X, t) 



fi(x,y,t) 



x=y 



+ -n(x, X, t)fi(x, X, t) 



(60) 



(To.) = . (61) 

As expected, the energy-momentum tensor matrix element is diagonal and can be ex- 
pressed in terms of the average energy density (e) and pressure (p). In four space-time 
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dimensions, we have 



{T^y) = diag ((e), a"^ (p) , a'^ (p) , (p) 



(62) 



Once more we can witness from ( |55| - pTD that infinities may appear in the hmit d = 3 because 
of the short distance behavior of the kernels Q and E. 

Now we discuss how to obtain the finite, renormahzed expectation value of T^i, in our 
Gaussian Ansatz, given by Eqs. (|ll]), (@), and (|49|) . First of all, we evaluate the dimension- 
ally regularized expression for (T^^), which is finite and well behaved. A nice feature of this 
regularization procedure is that it preserves the general covariance of (T^jy). Substituting 
(H) and (IID into (^^^), we obtain 

A . 1 



(Too) = - — TT^ + + ^ip^ + —a 
\ uu/ 2 2 ^ 4!^ 64 



d-lo2 T 

P -'2n+l/2 



+ 



1 



(47r)('^+i)/2(l -d) VAV 



A 



a 



'3-d 



2 2 2 I / 2 I A 2 \ 

YTdf +1^ +2^ J" 



2 (47r)('^+i)/2(l -d) VAV 



a y'^-3)/2^ 



(63) 



-2d 



+ 



-2 I ,„2 I A 4 X d-ln2T 

+ + 4!^ + ^ ^2n+l/2 



1 

64" 



(47r)(^+i)/2(l -c/) VA2; 



_ / 2 A 2\ _A 

r + 2*^ J " 2 (47r)('i+i)/2(l - d) va2; 



ci(l + c/) 

r 



1 



a 



(64) 



Next, we express the expectation value of the energy-momentum tensor in terms of the 
renormahzed parameters fiR and Xr. In order to do so, we substitute Eqs. ( ^6|) and 
into the last two expressions, resulting in 



n~° 1 



a 



(Too) = — vr^ + -a-^/3^J2„+i/2 + 



+ 



a 



647r2 



2A 



R 



2Xr 



(65) 
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„2 



647r2 



m 4 -i 



2Ar 2A 



(66) 



where we took the hmit d 3 and chose the scale A according to Eq. (^T]). 

One important feature of our Gaussian approximation is that the expectation value 
of Tfj_u turns out to be finite once we take into account the mass and coupling constant 
renormalization, and that n > 1. Therefore, the presence of interactions, more specifically 
the term ~g^u{X/8)QQ, leads to the cancelation of the infinities which exist in the free scalar 
field model fill . 



Now we discuss the limit Xpt = 0. In this limit, the renormalization prescription tells us 
that 

A = ; fil = fi\ (67) 

and hence we get back the unrenormalized free theory result for (T^z^) in which the diver- 
gences reappear from the 1/Xr terms. In order to have a well behaved free theory limit we 
must subtract the terms that diverge in the limit Aj? = 0. Moreover, this subtraction may 
be justified as a renormalization of coupling constants in a generalized Einstein equation 
provided that the entire subtraction is expressible in terms of covariantly conserved tensors 

In the free theory limit, the contribution — (7^,^(0; — /i|.)^/2AR vanishes since the structure 
equation of motion for a leads to a = fi^ for Aj? = 0. Therefore, we define the renormalized 
expectation value of the energy-momentum tensor as 

{T,.)n^{T,.)-g,.^. (68) 

Notice that this subtraction is basically a redefinition of the cosmological constant. 



V. CONCLUSIONS AND DISCUSSION 

In this work we used Dirac's variational principle and a simple Gaussian Ansatz to 
describe the time evolution of a self-interaction scalar field in flat Robertson- Walker space- 
time. Unlike ordinary perturbation theory, this approximation reflects some of the non-linear 
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features of the full quantum field theory which may be important for understanding various 
physical processes. Our trial wave functional was obtained from a general Gaussian state 
by choosing the momentum dependence of its kernels in such a way that we can recover the 
conventional vacuum in the limit that we have a Minkowski space-time. Nevertheless, we 
should point out that our vacuum state for a Robertson- Walker space-time is orthogonal to 
the one obtained in Refs. |p|,p!I|, where the Gaussian adiabatic vacuum was used, since they 
differ by terms of order 0(A;~^) for k oo. 

In principle, the equations for the variational parameters and the semi-classical Einstein 
equation may be used to study dynamical question about the universe, such as stability of 
the de Sitter space and the conditions for inflation setting in. In this case, however, we are 
confined to the chaotic inflation scenario since our Ansatz does not describe correctly the 
dynamics of the low momentum modes, as well as the Gaussian approximation in higher- 
dimensional field theory suffers from well-known shortcomings, analogous to the ones that 
appear in the the large- iV approximation. 
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